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PROPERTIES OF THE SCATTERING MATRIX AND 
DISPERSION ESTIMATES FOR JACOBI OPERATORS 

IRYNA EGOROVA, MARKUS HOLZLEITNER, AND GERALD TESCHL 


Abstract. We show that for a Jacobi operator with coefficients whose (j + 
l)’th moments are summable the J’th derivative of the scattering matrix is 
in the Wiener algebra of functions with summable Fourier coefficients. We 
use this result to improve the known dispersive estimates with integrable time 
decay for the time dependent Jacobi equation in the resonant case. 


1. Introduction 


This paper is concerned with Jacobi operators 

Hu(n) = a(n — 1 )u(n — 1) + b(n)u(n) + a(n)u(n + 1), 
satisfying a(n) -A 1/2, b(n) —> 0 such that 


a(n) 


+ \b(n)\e£l(Z) 




( 1 . 1 ) 

( 1 . 2 ) 


for some a > 1. Here £%(Z) is the set of all doubly infinite sequences for which the 
norm 


Me 


(E„ e z(i + HrK«)l p ) 1/p , pe[i,oo), 

sup„ ez (l + \n\) a \u(n)\, p = oo, 


is finite. The case a = 0 corresponds to the usual unweighted spaces £q = i v . 

The special case a(n ) = 1 is also known as the discrete Schrodinger operator. 
If (11.211 holds for a = 1 it is well known ([16]) that the spectrum of H consists of 
a purely absolutely continuous part covering [— 1,1] plus a finite number of eigen¬ 
values located in K. \ [—1,1]. In addition, there could be resonances at the edges 
of the continuous spectrum, which means that there exist corresponding bounded 
solutions (and which is equivalent to the fact, that the Wronskian of the two Jost 
solutions vanishes (cf. [5} Definition 3.5, Lemma 3.6])). In case of such a resonance 
it is a delicate question to determine the behavior of the scattering matrix near 
such a resonance under the minimal assumption a = 1. From the analogous result 
for the continuous Schrodinger equation it is expected that the scattering matrix is 
continuous at such a resonance and this was first established in [5] (see also [8] for 
a special case). In fact, in [3] this result was refined by showing that the scattering 
coefficients are in the Wiener algebra of functions on the unit circle with integrable 
Fourier coefficients. In the present paper we generalize this result by showing that 
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if G2* holds for a > 1 then also all derivatives of order up to a — 1 of the scattering 
coefficients are again in the Wiener algebra. 

This question is not only of interest in scattering theory but also plays an impor¬ 
tant role in the solution of the Toda equation via the inverse scattering transform 
(see e.g., muni where continuity of higher derivatives is needed) and in proving 
dispersive estimates for the corresponding linear evolution equation 

^ u(t) = Hu(t ). (1.3) 

The latter case has attracted considerable interest recently (e.g. HH3H3] and the 
references therein) due to its importance for deriving asymptotic stability of solitons 
for the associated nonlinear evolution equations (see e.g. .TtHT][T4j. 151). 

As an application of our results we will establish two dispersive decay estimates. 
First of all we extend [3J Theorem 5.3] to the case of Jacobi operators. 

Theorem 1.1. Let H be defined by (11.11) with \a — i| + |6| £ Then 

\\e~ ltH P c \\ e i_>e~ = 0{t~ 1/3 ), t ->• oo, (1.4) 

|| e — ltH P c \\ g 2 ^> ( 2 _" = 0(t -1 /2), t—> oo, a > 1/2. (1.5) 

The following extends [3] Theorem 6.1] to the case of Jacobi operators and is 
new in the case of resonances. 

Theorem 1.2. Let H be defined by (11.11) with |a — || + |6| £ K](Z), where j = 3 

if there is a resonance at z = +1 or z = — 1 and j = 2 if neither of these points 

is a resonance. Furthermore, for every resonance z £ {±1} let ipz{n) be a bounded 
solution of Hip zip) = — <Pz(n) normalized according to lim, w+00 (|y)f (n)| 2 + 

I Tz (— ri )| 2 ) = 2. Denote by P j the projection onto the span of <p~ given by the 
kernel [Pz](n, k) = ipz{ri)ipz(k) and set P~ = 0 if there is no resonance at z. By 
P ac we denote the projector on the absolutely continuous subspace of H. Then the 
following decay holds: 

1,6 ltHPac ~ = 0{t 4/3) ’ t ^ °°' (L6) 

Such a decay estimate with integrable time decay has previously only been es¬ 
tablished in the non-resonant case where P_i = P\ = 0 (cf. [3] Theorem 6.1]). For 
continuous one-dimensional Schrodinger equations an analogous result in the reso¬ 
nant case has been first proven by Goldberg [5]. For further results for continuous 
Schrodinger equations see el a- 

2. Properties of the Scattering Matrix 

In this section we look at scattering theory for the Jacobi operator H. As a 
general reference we refer to |T5J Chapter 10]. If (11.21) is contained in then 

there exist Jost solutions ip±(z,n ) of 

Hip±(z,n) = — <P±{z,n), 0 < \z\ < 1, (2.1) 

which satisfy \vca. n ^± 00 ip±( K z,n) = 1, where <p±(z,n) = ip±(z,n)z IfZn . These Jost 
solutions can equivalently be expressed as 

±oo 

ip±(z,n) = ^2,K±(n,£)z ±e -, 

i=n 


n £ Z, \z | < 1, 


( 2 . 2 ) 
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where the transformation operators satisfy 


\K±[n,i)\ < C±(n) 


±oo 

S(n,£) + (1-S(n,£)) Y 

fe=L^J 




(2.3) 


for ±£ > ±n and |_-J denotes the usual floor function. 

If furthermore ±n > =pl holds, we can replace C± (n) by some universal constant 
C±{n) < C. For 

±oo ±oo 

W±{z,n) = '^K±{n,l + n)z ±L = y^K±(n,£)z ±t 
1 =0 £=0 


we get similar estimates. Introduce the Wiener algebra 


-4 = {/(z) = Y f( m )z 


\\ e i < 00, \z\ = 1 


mE Z 


}■ 


with the norm ||/||_4 = ||/||^i. Formulas (12.21) and (12.31) imply 

<P±(z,n),(p±(z,n) £ A, (2.4) 

with || !p(z, u)||a independent of n for ±n > 0. Let us also recall the discrete version 
of the Wiener lemma [18J which states that if f(z) £ A and f(z) ^ 0 for all \z\ = 1 
then / _1 (z) € A. The estimate (12.31) immediately leads to the following property: 


Lemma 2.1. If (11.21) holds for a = j + 1 with j > 0, then -4-j(ip±(z,ri)) is an 
element of A for 0 < l < j. Moreover, for ±n > 0 the A-norms of these expressions 
do not depend on n. 


The fact that tp±(z 1 ,n) also solves the equation (12.11) gives rise to the scattering 
relations 


T(z)ip±(z,n) = R T (z)ip T {z,n) + ip T (z \n), \z\ = 1, (2.5) 

where the transmission coefficient T and the reflection coefficients R± can be ex¬ 
pressed in terms of Wronskians. To this end let us denote the discrete Wronskian 
by 

W ( f(z, n),g(z, n)) = a(n - 1) (/(*, n - 1 )g(z, n) - g(z, n - 1 )f(z, n )). (2.6) 

Introducing the functions 


W(z) = W{y+{z,l),V-{z,l)), W±{z) = W{ip T (z,l),(p±(z 1 ,1)), 
it follows that 

W±(z) 


T{z) = 


R±(z ) = =F- 


(2.7) 


2iW(z)’ ^ W{z) ' 

In [3j Theorem 4.1] it was proved that the transmission and reflection coefficients 
are elements of the Wiener algebra. Here we extend this result to derivatives. 


Theorem 2.2. If \a(n) — + |6(n)| £ £j +1 (Z), then fjfj(T(z)) £ A and fpjR±{z) £ 

A for 0 < l < j. 


Proof. We only focus on the resonant case W{\)W{— 1) = 0, since the other case 
is straightforward. Let z £ {±1} be a point with W(z) = 0. As a first step we 
introduce the expressions 

W±(z) = tp±(z,l)ip±(z,0) - <p±(z,0)ip±(z, 1). 
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Then the calculations in [5] Lemma 4.1] show 


W±(z) = C(*)*±(*), C (*) = —, (2-8) 

where 

zhoo 

*±(z)= J2 h±mzz) ± h ± (e) = ^\e)< f± (z,0)-<s> ( £\e)v±(z 1 i) ( 2 . 9 ) 

0— i±i 
2 

and 

±00 

$£ } (m) = ^ K±(k,£)z t . 

£=m 

Moreover, h±{t) satisfies a similar estimate as (12.41) . as we will show in Lemma [2741 
below. This immediately leads to the conclusion that ^±( 2 ) and all its derivatives 
up to order j are elements of the Wiener algebra. Next we need to distinguish 
between the cases tp + (z,0)ip-(z,0) 7 ^ 0 and tp + {z,l)ip-(z,l) 7 ^ 0 (these are the 
only cases since the solutions ip±(z,ri) cannot vanish at two consecutive points). 
Straightforward computations show 


W(z) 
z — z 





*-(*>- 


y-(^d) ,T, 

V+(z,l) ^ + 


(z 


V+(z,0)tp-{z,0) £ 0 , 

^+0M)<£-(M) 7^0. 


Since |T(z)| < 1, the zeroes of W(z) can at most be of first order, which shows that 
can vanish at most at —z. Hence if W (—5) 7 ^ 0 then € A by Wiener’s 
lemma (including all derivatives up to order j) and hence the same is true for T(z). 
Otherwise, if W(—z) = 0 and W(z) 7 ^ 0 one has the analogous properties for 
and using a smooth cut-off function (which is one near z and vanishes near —z) one 
can combine both results into one for and proceed is before. 

(z-z)(z+z) r 

To get similar results for R±, we use the following formulas: 


W ± {z) 
z — z 



<p+(z,0)<p-(z,0) ± 0 , 
V+{z,\)ip-{z, 1 ) / 0 , 


and proceed as in the previous case. 


□ 


Lemma 2.3. For h±(m) given by (12.91) we have 


±00 

|/i±(m)| < C J2 

n= im±lj 



+ \b(n. 


for some constant C > 0 and ±m > 0. 


Cr)±(m) 


( 2 . 10 ) 


Proof. We only consider the minus case here. First of all we need the discrete 
version of the Marchenko equation, which is given by ([11] Section 10.3]): 


±00 

I<± (■n , m) + ^2 K± i n > f)F± + m) 

£=n 


S(n , in) 
K±(n, n) ’ 


±m > ±n 


( 2 . 11 ) 






















SCATTERING AND DISPERSION ESTIMATES FOR JACOBI OPERATORS 


5 


where for F±(i) we have the following estimate: 


\p±m < c 


±oo 

E 

i=L|J 


( a{n) 



( 2 . 12 ) 


After some calculations, which can be found in [5J Lemma 4.1], the following equa¬ 
tion for h- can be obtained in the case <P-{z, l)y> + (£, 1) ^ 0 (the other case 
<P-(z, 1) = <p+(z, 1) = 0 is similar, cf. jj] Lemma 4.1]): 

h-{m)-J2 h -(t)v{t + m + l) = - <p-(z,l), (2.13) 

e=o "A ’ ' 

where v(£) = F_(£)z~ e . Now let —m > 2. We rewrite (12.1311 as 


h-(m) — ^ h-(£)v(m + £ + 1) = —H(m,N), where 

t=N 

H(m, N) = - 5°E f E -(+ !) + E h -( £ M e + m + 1) 

A_(°, 0) ^ 

and N < 0 such that < 1 with C given by (12.121) . The estimate 

\H(s,N)\ < C(N)rj-(m), m < -2, (2.14) 

follows from monotonicity of rj- and h_(-) £ £°°(Z_). Now we set 


/i_,o(m) = H(m,N), h_, k+1 (m) = ^ h- ik {£)v{£ + m + 1). 

t=N 


We show that 


|ft._ ifc (m)| < C(N)rj-(m) 


C E wm) , 


e=N 


(2.15) 


where C(N) is given by (12.141) . But this easily follows by induction, again using 
monotonicity of rj- and boundedness of /i_. □ 


For later use we note that in the resonant case the Jost solutions are dependent 
at z. If we define 7 via 

<P+(z,n) = jip-(z,n), (2.16) 

then a straightforward calculation using the scattering relations m as well as 
|T(fc )| 2 + |.R±(fc )| 2 = 1 shows 

T ^ = 7tE’ i2±(*)=±^4 (2-17) 

1 + 7 1 + 7 2 

In particular, all three quantities are real-valued since ip^{z,n) £ R and hence 
7 £ R. In the non-resonant case the scattering relations show 

T{z) = 0, R±(z) = - 1. (2.18) 

To establish Theorem ll.2l we also need the following generalization of Lemma l2.ll 

Lemma 2.4. Let (11.21) be contained in £j +1 (Z). Then g 

for 0 < l < j — 1. Moreover, for ±n < =pl, the A-norms of these expressions do 
not depend on n. 
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Proof. By (12.211 we have that 


<P±{z,n) - <p±(z,n) 
z — z 


±00 ^_|_£ ±oo ±£— 1 

E K±M Z = E E z k Z ±l ~ X ~ k 

e=±l £=±1 k—0 

±oo ±oo 

E( E 

e=o k=e± i 


Using this we conclude that the corresponding derivatives of the previous expres¬ 
sions have an .4-norm bound independent from n, if ±n < =pl by OK- □ 


In a similar way, we obtain the following lemma: 

Lemma 2.5. Let (11.21) be contained in £j +1 (Z). Then ag we n ag 

^■(^±( 2 )) are elements of A for 0 < l < j — 1, where ^±( 2 ) are defined in (12.91) . 

Proof. This follows as in the previous lemma using the estimate for h± from 
Lemma 12.31 □ 


Combining the last results we obtain: 

Theorem 2.6. Let (11.21) be contained in £j +1 (Z) and let 2 € {—1,1}- Then 

rf/rW-TWN £ ( *ro<i<i-i. 

az l \ 2 — 2 J dz L \ 2 — 2 J 


3. Dispersive Decay in the Discrete Case 

In this section we prove Theorem 11.11 and Theorem 11.21 For the one-parameter 
group of (0> the spectral theorem and Stone’s formula imply 


,-i tH 1 


27ri 


+ iO) — lZ(ui — iO)) duj, 


(3.1) 


where TZ(ui) = (H — u>) 1 is the resolvent of the Jacobi operator H and the limit 
is understood in the strong sense H3 Problem 4.3]. For the kernel of the resolvent 
R(z) = {H — —) _1 , we have the following explicit formula for 0 < | 2 | < 1 (cf. 

[Ill (1-99)]): 


[R{z)\{n, k ) = 


1 [ ip+(z,n)tp-(z 1 k) for n>k, 


W(z ) 1 tp + (z,k)T-(z,n) for n<k. 


(3.2) 


Formulas m and (13.21) then lead to the following explicit representation of the 
kernel of the propagator e~' tH P ac - 


[e~ itH Pac\ ( n,k ) = ^f e- itcose v + (e w ,n)<p-(e ie ,k)T(e ie )d9. (3.3) 


We also need a small variant of the van der Corput lemma, which can be found in 
[31 Lemma 5.1]: 
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Lemma 3.1. Consider the oscillatory integral I(t) = e ltv ^ f(9)dd with [a,b\ C 
R some compact interval and v(6) real-valued. Let min |i/ J )(d)| = mj > 0 for 


some j > 2 and let 


f(0) = /( m ) e ‘ m<? f or 6 e \ a M with ll/lli = 5Z I/Ml < 00 • ( 3 - 4 ) 

m(=Z m(zL 


Then 


IMI < 


Ml/lli 

{mjtf/j ’ 


where Cj is an universal constant. 


fort > 1 , 


(3.5) 


The proofs of Theorem ll.ll and the non-resonant case of Theorem 1 1.2 1 can now be 
done in exactly the same way as in [3] Theorem 5.3, Theorem 6.1]. Thus it remains 
to consider the resonant case of our main Theorem o As a first step we have a 
closer look at our projection operators in the next lemma: 


Lemma 3.2. The following expressions for the kernels of our Projectors P~ are 
valid: 




[Pi]{n,k) = <p+(l,n)p_(l,fe)T(l) 


1 C 2 j 

I e -itcos e de __ 


27T 


\/27rit 


= [P_i](n,fc) = <£+(—1, ro)</5_(—1, k)T(—l) 


i 


itn 


e~ itcos9 de 


+ 0(t~i), 
(3.6) 

1 n 


it7T 


+ Q(t~ ,. 


Proof. It is clear that ip£ (n) = c±tp± (z, n ) and by our normalization c_ 2 + c + 2 = 2. 

Using (12.161) we have c_ = 70 + and hence c± = 1 + (J T . Moreover, (12.171) implies 
c + C-T(z ) = 1 and hence [P~\(n,k) = ip + (z,n)ip-(z,k)T(z). Furthermore, by [12l 
10.9] and [T2] 11.5], we get 

i(J 0 (t)-i^o(t)) = 2 - [[e- itcose d6, 

where Jo(t) denotes the Bessel function, and Ho(t) the Struve function, both of 
order 0 (cf. [12] 10 . 2 ] and |T2] 11-2]). Moreover, by [12; 11.6] the following asymp¬ 
totics hold 

Ho{t)-Y 0 (t) = — +O(t- 3 ), 

7 Tt 

where Yo(t) denotes the Neumann function of order 0 (cf. [12] 10.2]). Since by 
[12] 10.4] H^\t) = Jo(t) + i Yo(t) and HQ 2 \t) = J 0 (f) — iTb(i), the claim follows 
using the asymptotics [12j 10.17.5, 10.17.6]. The argument for the second formula 
is analogous. □ 


Since we have got everything we need,we can finish the proof of our main theorem 
now: 


Proof of Theorem \1.2\ Using (13.31) . we get 

[e- itH P ac ](n,k) = 7 ^ f e~ itcose G(e,n,k)de , 


tojco 
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where 

G(0,n,k ) = e i|n - fe| V+(e iS ,fc)^-(e ie ,n)T(e ie ). (3.7) 

Taking advantage of the fact that all functions here are 27r-periodic, we can integrate 
along the interval [—f, %-] instead. By Lemma [3721 


V(n,k,t) := [e ,tH P a c - -j==rjPi - P-i](n,k) 




it cos 6 


( G(9 , n , k ) — G(0, n, k))dd + G(0, n, fc) 7 


itn 


(3.8) 




i t cos 6 


(G(9 , n, k ) — G(7 t, n, k))dd — G(7r,n, k )-1- 0(t *). 

ltj[ 


(3.9) 


If we integrate (13.81) by parts, we obtain 


1 

27T 



3~ itcose (G(e, n ,k) - G{0, n ,k))d0 + G(0,n,k)J- = -Vg 

x ' ltn 2mt 






1 [* e -itcose_^_ / G(6>, n, k) - G(0,n,fc) 

27rif J_* \ sind 


If we do the same in (13.91) and use G(^-,n,k) = G(—^,n,k) we see that all the 
terms of order j vanish and hence 


V(n, k, t ) 


i f 2 —it cos o d 

2nt J d6 

i j 2 —it cos 9 d 

2nt Jn d6 


( G(0, n, k) — G(0, n, k) \ ^ 
\ sind J 


(3.10) 


G(9, n, k ) — G(n, n, k) 
sin 9 


d6 + 0{t~i), 


where we neglect the summand of order t 2 from now on. The other two summands 

4 

are treated separately and we show the desired t~ 3 time decay for each expression. 
Since the calculations are similar, we only focus on (13.101) . Next, to apply Lemma 
I3T1 we split the domain of integration into parts where either the second or third 
derivative of the phase —it cos 6 is nonzero. This gives us the time decay. It remains 
to show that the derivatives with respect to 6 are elements of A with „4-norms at 
most proportional to (|n| + |fc|). Of course this will in general not be true because 
of the possible singularity at the other end of the continuous spectrum. However, 
since this possible singularity is outside the domain of integration we can redefine 
our functions there. 

We distinguish the cases (i) n < 0 < k, (ii) 0 < n < k, and (iii) n < k < 0. 
Introduce the functions 


m = 


e ie - 1 
sin 9 ’ 


g(9,n, k) = G(0, n, fc) e - i|ri " fc|e , 


h(9, n, k) 


g(0, n,k) - g(0,n,k) 
e ie — 1 


Then 

d G(9,n,k) — G(0,n, fc) 
dO sin 9 


w I ^ 


gi (k-n)9 _ 1 


cXO 


- 1 


g(6,n,k) + f(6)h(9,n,k) 
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Since only the values of / for 8 £ [— ^are relevant, we redefine it on [^, 
such that f, f £ A (e.g. by multiplying with a smooth periodic function which is 1 
on [—§,§] and vanishes near 7 r). 

We observe that 


gi (k—n)6 _ -y gk—n _ y 


k—n—1 


pi 9 


-1 


z-l 


= E z " eA 


t =0 


with its .4-norm bounded by k — n, and that of its derivative by (k — n ) 2 . Now it re¬ 
mains to consider the M-norms of g(9, n, k ), j§g(9, n, k ), h(9 , n, k), and j§h(9, n, k). 

We start the with case (i). Then g(9,n,k) G -4 with M-norrn independent of 
n and k. After applying the product rule, Lemma 12.11 and Theorem 12.21 imply 
\\-^9(9,n,k)\\ A < C. Moreover, by (ESl 

h(9, n, k) = --- l P+{z, k)<p-(z, n) +- - - ip-(z,n)T( 1) 


+ 


z — 1 
<P-(z,n) - <p- (1, n) 
z-l 


<P+(z,k)T(l), 


z-l 

je 


z = e 


Invoking Lemma, 12.11 Lemma [2~4l and Theorem 12.61 proves that h and its deriv¬ 
ative are again contained in A, with ,4-norms independent of n and k , thus we are 
done in case (i). In the cases (ii) and (iii) we use the scattering relations <123)1 to 
get the following representations: 


g{9,n,k) = 


<P+(z,k) (R+(z)fi + (z,n)z 2n + ip+(z \n)) , 0 < n < k, 
<P-(z, n) ( R-(z)ip-(z, k)z 2k + k)) , n < k < 0, 


where z = e' e . Thus g(9,n,k) has an M-norm independent of n and k, since for 
any f € A and to £ Z, we have that f(z)z m £ A , with M-norm independent of to. 
Taking derivatives with respect to 9 we get additional terms from the derivatives 
of z 2n and z 2k , so by Lemma [2~T1 and Theorem 12.21 it follows that \\-^g{9,n,k)\\A 
is at most proportional to |n| and |/c| respectively. Finally it remains to consider 
h(9,n,k). In case (ii) h(9,n,k) can be represented as follows: 


9.(9, n, k ) = 


<P+(z,k) - £+(!,&) 


z-l 
R + (z)-R + (l) 


tp+(z 


-l 


(p+(z ■ ,n) - <p+(l,n) 


z-l 


<p+(l,k) 


k 


+ 


Z — 1 

<P+(z,n) - p+(l,ra) 


tp + (z,n)(p + (z,k)z 2k ■ 
(p+(z,k)R+(l) ■ 


- 1 


R+(l)<P+(z,n)(p+(z,k) 

<P+{l>n)R+(l), 


z — 1 
<P+(z,k) - &+(l,k) 


z-l 


for 2 = - ie 


z-l 

e‘“. Here again every summand is an element of A by Lemma 12.11 
Lemma m and Lemma 12.61 Since the derivative of p also occurs here, 
we conclude that the ,4-norm of -^h(9,n,k) is at most proportional to \k\ 2 . In 
the case (iii) this derivative will be proportional to |n| 2 . This indeed proves the 
claim about G ( 9 ’ n ’ k )- G (°’ n < k ) f or 9 £ [— f j §]• Since the same procedure works for 
G(8,n,k)-G(ir,n,k) ^ q g [§, ^-] this finishes the proof. □ 
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